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Abstract 

Introduction: In the departments of the Vaupes and 
Guaviare, in southeastern Colombia, in a transitional 
area between Amazonia and the eastern plains, inhabit 
indigenous groups belonging to the Tukanoan (East) and 
Guahiban linguistic families. Although some studies have 
dealt with the culture and the cosmology description of 
these groups, little research has been done on the biological 
diversity and genetic relationships of such groups. 
Objective: To estimate the diversity, the structure, and the 
genetic relationships of one Guahiban and two Tukanoan 
groups of the Colombian Amazonian region. 
Methods: Samples were collected (n= 106) from unrelated 
individuals belonging to the Vaupes native indigenous 
communities. The DNA was extracted and nine autosomal 
microsatellites were typed. Several measures of diversity, 
FST, pairwise FST, and population differentiation between 
groups were calculated. Finally, it was estimated the genetic 
distances of the groups studied in relation with other 
Amazonian, Andean and Central American indigenous 
people. 

Results: 1. The genetic diversity found stands within 
the range of other Amazonian populations, whereas 
compared to the mestizo and afro-descendant Colombian 
populations, such diversity showed to be lower. 2. The 
structure and population differentiation tests showed 
two clusters; one consisting of the Vaupes Tukanoan and 
Guaviare Tukanoan groups, and a second one formed by 
the Guayabero. 3. Tukanoan groups are found to be closer 
related to the Brazilian Amazonian populations than to the 
Guayabero. 

Conclusion: The results of this study suggest that 
the Guayabero group from Guaviare, are genetically 
differentiated from those Tukanoan groups of the Vaupes 
and Guaviare. 



Resumen 

Introduccion: En los departamentos del Vaupes y el 
Guaviare, al Suroriente de Colombia, en un area transicional 
entre la Amazonia y los Llanos Orientales, habitan pueblos 
indigenas de las familias linguisticas Tucano (Oriental) 
y Guahibo. Se han hecho estudios sobre la cultura y la 
cosmologia de dichos grupos, pero poco se ha investigado 
sobre su diversidad biologica y sus relaciones geneticas. 
Objetivo: Estimar la diversidad, la estructura y las relaciones 
geneticas de dos poblaciones Tucano Oriental y un grupo 
Guayabero (Guahibo) de la Amazonia colombiana. 
Metodos: Se recolectaron muestras de 106 individuos no 
emparentados, pertenecientes a comunidades indigenas 
del Vaupes y el Guaviare. Se tipificaron nueve microsatelites 
autosomicos. Se estimo la heterocigosidad, la diversidad 
genica, FST y FST pareados entre grupos. Se realizo una 
prueba de estructura espacial y se construyo un arbol 
de distancias geneticas con otras poblaciones indigenas 
amazonicas, Andinas y centroamericanas. 
Resultados: La diversidad genetica se encuentra en el rango 
estandar de otras poblaciones amazonicas, pero es mas baja 
que en poblacion mestiza y afrodescendiente de Colombia. 
Los Tucano Oriental del Vaupes y el Guaviare conforman 
una unidad genetica que se diferencia significativamente 
de los Guayaberos. Los grupos Tucano estan mas cerca 
geneticamente de las poblaciones amazonicas brasileras 
que de los Guayabero. 

Conclusion: Los resultados de este estudio sugieren 
que los Guayaberos en el Guaviare, estan diferenciados 
geneticamente de los grupos Tucano del Vaupes y de los 
Tucano del Refugio en el Guaviare. 
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Introduction 

The Vaupes and Guaviare regions in southeastern 
Colombia are inhabited by several indigenous groups 
belonging to different linguistic families. These groups 
have an extensive history of demographic expansions, 
wars, alliances, migrations and trades. With the arrival of 
European populations into the region, a new chapter was 
opened, leading to the disruption of some cultural and 
social dynamics of the existing indigenous groups. To a 
large extent, it can be said that such disruptions caused a 
major reduction on the size of the indigenous populations, 
displacements from original territories, or even a total 
extinction of some groups 1 . 

Currently, 24 ethnic groups live in this area representing 
about 30% of the national linguistic diversity. These groups 
have been classified 2 into five linguistic families: Tukanoan 
(represented by the Eastern Tukanoan subfamily), 
Guahiban, Arawakan, Kakua-Nukak and Cariban. These 
ethnic groups are found scattered throughout both the 
Vaupes and Guaviare departments, East Tukanoan groups 
being the most widely distributed in the Vaupes, while in the 
Guaviare highlights predominates the presence of Guahibos, 
with some presence of Eastern Tukanoans whom in recent 
years have migrated from the Vaupes. Most of these ethnic 
groups live in indigenous communities with populations 
ranging from 22 to 6,222 inhabitants 3 and some of them are 
currently in serious danger of disappearing, as is the case of 
the Eastern Tukanoan Pisamira in the Vaupes 2 . 

Great efforts have been made towards an understanding 
of the social and cultural dynamics of the groups from 
this region 4 7 . However, very little (and very fragmented) 
information have been obtained at the genetic level 8,9 . 
With the aim to expand our biological knowledge on the 
populations of the Colombian Amazonia, nine autosomal 
microsatellites were typed to evaluate diversity, gene flow, 
genetic structure and degree of genetic differentiation 
between Eastern Tukanoan and Guahibo (Guayabero) 
native groups of the Vaupes and the Guaviare. 

Materials and Methods 

Samples. For this study were collected biological samples 
of peripheral blood (4 mL) of 106 individuals belonging to 
indigenous groups in the departments of the Guaviare and the 
Vaupes, in southeastern Colombia (Fig. 1). The participants 
were healthy adults (men and women) without biological 
relationship (until the third degree of consanguinity). 
Samples were obtained with previous informed consent from 
each participant and with the agreement of the community 
leaders. This study was approved by the Institutional Human 
Ethics Review Committee at the Universidad del Valle, 
following the guidelines of the resolution 008430 of 1993 of 
the Ministry of Health of the Republic of Colombia. 
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Figure 1. Map of continental Colombia with the location of the 
sampling sites: MM (Tucanos from Vaupes), El Refugio (Tucanos 
from Guaviare), and el Barrancon (Guayaberos from Guaviare) 




Populations. The sample from Mini (n= 55), named here 
Vaupes Tukanoans, groups Tukanoan 2 speakers, some of 
whom live in the municipality of Mitu(ri5'29" N, 70°14'10" 
O) and others are more like a "floating population" arriving 
in Mitu to trade goods and food but whose villages are 
elsewhere. The sample of El Refugio (2°34'07" N, 72°35'54" 
W) in Guaviare (n= 22), named here Guaviare Tukanoans, 
groups Tukanoan speakers who migrated from the Vaupes 
about 50 years ago. The sample of El Barrancon (2°15'29"N, 
70°14T0" W) in the Guaviare (n= 29), named here Guaviare 
Guayaberos, groups people speaking Guayabero, classified 
within the Guahibo 2 language family. 

Laboratory. The DNA was extracted by the salting out 
procedure 10 . Each microsatellite (TH01, D7S820, D13S317, 
vWA, FGA, TPOX, FOB, LPL, FES / FPS) was amplified by 
PCR using primers whose sequences were obtained from 
STRBase 11 . The fragments were separated by polyacrylamide 
gel electrophoresis (8%), stained with silver nitrate. Each 
individual genotype was determined using allelic ladders 
constructed from DNA of previously established by capillary 
electrophoresis. Data analysis. The allele frequencies 
were obtained by counting; the intralocus independence 
(Hardy- Weinberg equilibrium) was tested. The observed 
heterozygosity (Ho) of each locus was calculated as the 
proportion of heterozygotes in each population. The 
observed heterozygosity of each population is the average 
heterozygosity across all loci. The gene diversity for each 
locus was cal culated as (n / (n-1)) (1-Sp. 2 ), where n is the 
sample size and p. is the frequency of allele i. 
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The total gene diversity was calculated as the average 
diversity over all loci (assuming panmixia, genetic diversity 
is equivalent to the expected heterozygosity He). The 
calculations were made with the software Arlequin 3. II 12 . 

Through an analysis of molecular variance, the degree 
of population genetic structure in two hierarchical 
levels was estimated: among groups (Vaupes Tukanoans, 
Guaviare Tukanoans and Guaviare Guayaberos) and 
among individuals within groups. This analysis divides 
the total variance into covariance components due to the 
intraindividual, interindividual and/or interpopulation 
differences. Covariance components were used to compute 
fixation indices, F . Gene flow was measured as Nem = (1/ 
F sr -l)/4. The calculations were performed with Arlequin 
3.11. 

To determine the degree of differentiation among the three 
populations studied, exact test was performed based on the 
distribution of the genotypic frequencies, the calculations 
were made in GENEPOP version Web 13 . The Geneland 14 
software was used to implement a model of spatial and 
correlated allele frequencies and the number of populations 
K was estimated (with a maximum number of populations 
of 10 and 100,000 iterations), this software take into account 
the genotypic data and geographical position of the samples, 
regardless of linguistic affiliation. 

The allelic frequencies of four loci (D7S820, D13S317, TH01, 
TPOX) were compared with data sets from six Amazonian 
populations, seven Andean and two Panamanian 
populations previously reported 1518 . Genetic distances D A 
were calculated and graphically depicted with a Neighbor- 
joining tree (NJ) using the program POPTREE2 19 . The 
consistency was tested by bootstrap. 

Results 

Allelic frequencies. The distribution of allelic frequencies 
of the Vaupes Tukanoans, Guaviare Tukanoans and 
Guaviare Guayaberos are shown in Table 1. The FGA was 
the microsatellite with the largest number of alleles and 
higher values of observed heterozygosity in the Guaviare 
Tukanoans and Guaviare Guayaberos samples (Table 2). The 
LPL was the system with the lowest number of alleles and 
observed heterozygosity in the Vaupes Tukanoans and the 
Guaviare Tukanoans. The FGA exhibited a deviation from 
the Hardy- Weinberg equilibrium in the Vaupes Tukanoans 
and Guaviare Guyaberos while systems D13S317, vWA, 
FOB, showed a significant deviation from HW equilibrium 
in Guaviare Tukanoans. Private alleles (relative to the three 
groups) were obser ved in the Vaupes Tukanoans at loci 
THOl (allele 8), D13S317 (15), FGA (20), TPOX (9) and in 
the Guaviare Guayaberos at loci TPOX (7), FOB (6), FES/ 
FPS (8). All these alleles have been previously reported for 



Table 1. Allelic frequencies in nine autosomal STRs in three Amerindian popula- 
tions of the Vaupes and the Guaviare departments, Southeast Colombia. 



Locus and alleles 


Vaupes 
Tukanoans 


Guaviare 
Tukanoans 


Guaviare 
Guavaberos 


THOl 








6 


0.323 


0.361 


0.478 


7 


0.49 


0.417 


0.522 


8 


0.021 


0 


0 


9.3 


0.166 


0.222 


0 


Chromosomes 


96 


36 


46 


D7S820 








9 


0.01 


0.028 


0.048 


10 


0.156 


0.417 


0.262 


11 


0.438 


0.222 


0.214 


12 


0.323 


0.277 


0.428 


13 


0.073 


0.056 


0.048 


Chromosomes 


96 


36 


42 


D13S317 








9 


0.312 


0.111 


0.75 


10 


0.219 


0.25 


0 


11 


0.146 


0.056 


0 


12 


0.063 


0.166 


0.1 


13 


0.104 


0.056 


0.125 


14 


0.146 


0.361 


0.025 


15 


0.01 


0 


0 


Chromosomes 


96 


36 


40 


vWA 








14 


0.01 


0.083 


0.025 


15 


0.115 


0.111 


0 


16 


0.469 


0.416 


0.85 


17 


0.271 


0.278 


0.1 


18 


0.104 


0.056 


0.025 


19 


0.031 


0.056 


0 


Chromosomes 


96 


36 


40 


FGA 








19 


0.115 


0.063 


0.025 


20 


0.01 


0 


0 


21 


0.021 


0.094 


0 


22 


0.042 


0.031 


0.075 


23 


0.094 


0.031 


0.125 


24 


0.25 


0.25 


0.075 


25 


0.208 


0.187 


0.35 


26 


0.208 


0.25 


0.15 


27 


0.042 


0.094 


0.175 


28 


0.01 


0 


0.025 


Chromosomes 


96 


32 


40 


TPOX 








7 


0 


0 


0.023 


8 


0.404 


0.305 


0.727 


9 


0.032 


0 


0 


10 


0.021 


0.056 


0 


11 


0.319 


0.417 


0.25 


12 


0.224 


0.222 


0 


Chromosomes 


94 


36 


44 


F13B 








6 


0 


0 


0.022 


8 


0.011 


0 


0.022 


9 


0.076 


0.083 


0.109 


10 


0.38 


0.333 


0.304 


11 


0.511 


0.556 


0.5 


12 


0.022 


0.028 


0.043 


Chromosomes 


92 


36 


46 


LPL 








10 


0.848 


0.778 


0.295 


11 


0 


0.028 


0.046 


12 


0.152 


0.194 


0.659 


Chromosomes 


92 


36 


44 


FES/FPS 








8 


0 


0 


0.031 


10 


0.089 


0.088 


0.063 


11 


0.622 


0.706 


0.406 


12 


0.167 


0.177 


0.406 


13 


0.089 


0 


0.094 


14 


0.033 


0.029 


0 


Chrntnnsntnes 


90 


34 


3? 
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Table 2. Genetic diversity indexes of three Amerindian populations of the Vaupes and the Guaviare. 







Observed 




Average number of 


Populations 


n 


heterozygosity [H n ) 


Gene Diversity 


alleles 


Vaupes Tukanoans 


48 


0.665 +/" 0.199 


0.622 +/" 0.333 


5.444 +/" 2.186 


Guaviare Tukanoans 


18 


0.663 +/" 0-205 


0.632 +/" 0.345 


4.778 +/" 1-641 


Guaviare Guayaberos 


23 


0.544 +/" 0.179 


0.466 +/- 0.266 


4.444 +/- 1.810 



other Andean and Amazonian populations. Two individuals 
of the Vaupes Tukanoans presented TH01 allele 8, which 
had not been reported for this group. 

Genetic diversity. The observed heterozygosity of Tukanoans 
and Guayaberos was found in the range 0.665-0.544, the 
gene diversity in the range 0.622-0.466 and the average 
number of alleles in the range 5,444-4,444 (Table 2). The 
Guaviare Guayaberos presented the lowest values of genetic 
diversity, although not significantly different from those 
found in Tukanoan groups. 

Genetic structure and population differentiation. The F 
estimated by AM OVA for the three groups was 0.06691 
(p-value <0.00001). The F estimated across all loci was 
0.07477 (IC 95% 0.02897-0.13481). With respect to the 
pairwise F ST (Table 3), the Guayaberos showed moderate 
differentiation with high significance (p <0.01) regarding 
the two Tukanoans groups, while the FST between the two 
Tukanoan populations showed no significant differentiation 
(p> 0.05). 

The exact test of population differentiation showed highly 
significant differences between Guayaberos and the two 
Tukanoans groups (p <0.01). However, it was not possible 
to detect any significant differences between the Vaupes 
Tukanoans and Guaviare Tukanoans, showing that genotypic 
frequencies of each locus havea similar distribution in the 
two populations. Similarly, the test of spatial differentiation 
in Geneland showed in all cases (ten separate runs to test 
the consistency of the results), the separation between 
Guayaberos and Eastern Tukanoan groups, with the two 
Tukanoan populations forming a single undifferentiated 
population. 

N m values show a low number of migrants between 
Guayaberos and the two populations of Eastern Tukanoan 
(Table 3), the opposite happens between the two Tukanoan 
populations showing a high rate of migration. 

Andes- Amazonia Genetic distances. The NJ tree (Fig. 2) 
shows two clusters East-West: the Amazonian Tukanoan 



populations on a branch and Guayaberos, Andean and 
Panamanian groups on the other. 

Discussion 

The genetic diversity of the Tukanoan populations (0.62 to 
0.63) stands within the range reported for other Amerindian 
South American populations 16 (0.554-0.753), which is 
lower than the african descendant or mestizo Colombian 
populations 15,20 . These low diversity values compared to those 
of European and African human groups havebeen extensively 
reported 21 for other autosomal markers, Y chromosome, 
and mitochondrial DNA. Such low values could be related, 
first, with the peopling of America (by small populations 
that spanned Bering Strait to North America about 15,000- 
25,000 years ago, reaching South America around 15,000 
years ago 22 , and second, with historical and demographic 
events that have caused the redution in population sizes since 
the arrival of the Spaniards. The Guayaberos presented the 
lowest diversity indices, which although is not segnificantly 
different from the indices Tukanoans, might reflect a different 
biological history between these two groups. The presence of 
four private alleles in the Vaupes Tukanoans and three alleles 
in Guayabero can be interpreted as an additional evidence 
of the separation between these two groups. The lowest 
diversity of Guayaberos, with an expected heterozygosity 
(0.466) outside of the range found in other South American 
Amerindian populations (0.554-0.753), would be associated 
with historical events typical of nomadic and seminomadic 
groups, coupled with the loss of alleles causedby the reduction 
of their population sizes as a consequence of wars with other 
groups in the region in prehispanic times, in addition, some 
possible infectious diseases carried by the non-native people 
caused a massive amount of deaths among the Amerindian 
populations during the first half of last century 23 . Also, more 
recently, the pressure of illegal armed groups have led to 
not only death of many indigenous peoples but also to the 
territorial displacement of such groups from their villages on 
the Ariari and Guayabero rivers to the current population 
of Barrancon 24 , destroying much of the cultural and genetic 
structure of these communities. 

Additionally to the reduction of population sizes, another 
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Table 3. Pairwise FST values of the three Amerindian populations of the Vaupes and the Guaviare. 



Pair of Population 




p 

r ST 


p Value 


N m 


Vaupes Tukanoans - 


Guaviare Tukanoans 


0.00752 


0.07207±-0.0227 


32.99 


Vaupes Tukanoans - 


Guaviare Guayaberos 


0.09098 


0.00000±0.00000 


2.5 


Guaviare Tukanoans 


- Guaviare Guayaberos 


0.12341 


0.00000±0.00000 


1.78 



important factor in the structuring of the genetic variation in 
Guayaberos is the inbreeding marriage practices 25 . For these 
groups, kinship represents the major element of identity, 
and the marriages occur between bilateral cousins, a factor 
that contributes to the loss of genetic diversity The Vaupes 
Tukanoans by contrast, associate to practices of linguistic 7 
exogamy, in which marriages are to occur between speakers 
of different language groups. In this exogamic maritalnet, 
women migrate to the community of their future husbands, 
allowing a constant gene flow throughout the region. It is 
possible that this practice has helped to maintain average 
levels of diversity in populations with sizes not exceeding 
one thousand people, as in the vast majority of ethnic 
groups in the Vaupes 3 . However, a more detailed study with 
uniparental markers is needed to allow for a more in depth 
evaluation of the impact that matrilocality/patrilocality 
oriented practices, asymmetry of gene flow, and inbreeding 
have had on the biological structure of the ethnic groups in 
the Vaupes region. 

The Exact Test of population differentiation, the spatial 
structure test, and pairwise F sr , support the biological 
separation between the Guayaberos and Eastern Tukanoan 
groups. Culturally, Guayaberos have not participated, nor 
they currently do, in the exogamous system of the Vaupes. 
This lack of exogamous marital net is reflected in the results 
presented here. On the other hand, geographical proximity 
is not synonymous of a common biological history The 
Guaviare Eastern Tukanoans are located near the Guayaberos 
(~ 7 km) ever since approximately five decades ago. 

Nonetheless no evidence of marriages or a high gene 
flow between these communities were reported during 
the field work period carry out for this study The 
estimated number of migrants between the Tukanoan and 
Guayabero groups was low compared to the number of 
migrants among the Vaupes Tukanoans and the Guaviare 
Tukanoans. This can potentially be a good evidence 
that in human populations socio-linguistic practices are 
crucial in the structuring of genetic variation; this can 
contribute or limit the gene flow between populations. 

Comparison with other indigenous populations through 
the estimation of genetic distances, group Tukanoans in the 
cluster of Amazonian populations. This grouping should 



not be interpreted as a common origin; nonetheless it 
could represent an ancient relationship among Tukanoans 
groups and other Brazilian Amazonian populations; these 
relationships are reflected in both oral tradition and cultural 
practices. The Tukanoan groups, for example, place their 
mythological origin in the "milk lake", an unspecified 
location downstream the Amazon River 6 . From there, an 
ancestral anaconda transported them (Tukanoans) upstream 
along the Vaupes river, and distributed the groups along 
the route. Such mythological story of origin is consistent 
with the genetic data analyzed here, which would support 
a possible origin of Tukanoan groups further southeast in 
the region now occupied by these groups in the Brazilian- 
Colombian border. 

The case of Guayabero and Andean populations is similar to 
Tukanoan and the Amazonian populations. It is not possible 
to trace with certainty an origin of the Guayaberos, but the 
evidence here presented suggests that these populations 
may possibly have an Andean-Eastern plains origin more 
rather than an Amazon origin, notwithstanding some of 
the Guayabero's ancient practices — such as consumption 
of manioc — which presents some resemblances to the 
Tukanoan and Arawakan practices, in the Amazon Basin. 

The Guayabero location in ancient times (on the headwater 
of the Ariari and Guayabero rivers 25 , in the foothills of the 
Cordillera Oriental) and his nomadic condition, would 
have allowed them a greater mobility along the foothills 
and forests of the Orinoco savanna, thus presumably having 
more contact with Andean indigenous than with Amazonian 
groups. 

Despite the genetic differentiation found between 
the Guayabero and Eastern Tukanoans groups, the 
interpretation of the East-West genetic relationships of 
these groups should take into account the limitations that 
carry the type and number of loci examined, sample sizes, 
and number of populations involved here. Interdisciplinary 
studies providing ethnographic, linguistic, and biological 
information of a higher number of Tukano, Guahibo, 
Arawak, Kakua-Nukak, and Puinave populations of the 
Eastern Plains and the Amazonia, are needed in order to 
establish a more reliable tree of the genetic relationships of 



137 



Braga Y. et al /Colombia Medica - Vol. 43 N°2, 2012 (Apr-Jun) 



Figure 2. Neighbor-Joining tree based on distances (DA) of 16 South American Amerindian populations and two Panamanian populations (Embera, 
Ngobe). The numbers in the nodes indicate the bootstrap value. Current study: Tucanos- Vaupes, Tucanos-Guaviare, Guayaberos. 



60 



13 



13 



I — Vaupes Tucanoans 
— Guaviare Tucanoans 



Zoe 



63 



18 



40 



Awa-Guaja 



Urubu-Kaapor 
- Parakana 



Wayampi 



1 — Tiriyo 

72 



East - Amazonia 



28 



18 



20 



Guayaberos 



1 — Ngobe 

- Matsiguenga 

- Cayapa 



22 



i — Quechua 



27 



Quechua-Matsiguenga 
Pijao 



55 



52 



92 



Embera 



Yanacona 



Paez 



West - Andean 



0.02 



these groups. 
Conclusions 

1. The genetic diversity of Tukanoans is within the range of 
values of other Amazonian indigenous populations, which 
is lower than mestizo and Afro-descendant Colombian 
populations. 

2. The test of population differentiation, the pairwise F ST 
138 



indices, and the test of spatial structure in geneland do not 
detect significant differences between the Vaupes Tukanoans 
and Guaviare Tukanoans. 

3. The same tests detected significant differences between 
Eastern Tukanoans and Guayaberos, so it is concluded that 
although these two groups belonging to the Amazon region 
or transition to the eastern plains, have had an isolation level 
such today they constitute differentiated populations. 

4. While there is a short distance between the community 
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of El Refugio (Guaviare Tukanoans) and the community of 
El Barrancon (Guaviare Guayaberos), culture and language 
would be acting as a barrier to gene flow, maintaining a level 
of genetic differentiation between these two groups. 
5. The Tukanoans are located closer to the populations of the 
Brazilian Amazon than the Andean populations. Opposite 
happens with the Guayabero, showing more proximity to 
the Andean populations than the Amazon. 
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